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ABSTRACT 



Slim-disk models describe accretion flows at high luminosities, while reducing to the standard thin disk form in the low luminosity 
limit. We have developed a new spectral model, slimbb, within the framework of XSPEC, which describes fully relativistic slim-disk 
accretion and includes photon ray-tracing that starts from the disk photosphere, rather than the equatorial plane. We demonstrate the 
features of this model by applying it to RXTE spectra of the persistent black-hole X-ray binary LMC X-3 . LMC X-3 has the virtues 
of exhibiting large intensity variations while maintaining itself in soft spectral states which are well described using accretion-disk 
models, making it an ideal candidate to test the aptness of slimbb. Our results demonstrate consistency between the low-luminosity 
(thin-disk) and high luminosity (slim-disk) regimes. We also show that X-ray continuum-fitting in the high accretion rate regime can 
powerfully test black-hole accretion disk models. 

Key words, black holes - accretion disks - LMC X-3 



1. Introduction 

Accretion powers the most luminous objects in the universe by 
converting gravitational potential energy into highly energetic 
radiation. In order to understand this mechanism it is crucial to 
study the primary engine of accreting black-hole systems: the 
accretion disk. Accretion disks form when gaseous matter, usu- 
ally an assembly of free electrons and various types of ions, spi- 
rals onto a central gravitating body by gradually losing its initial 
angular momentum as a result of viscous and magnetic stresses. 
The simplest ana lytical model of an accretion disk is the standard 
thin d isk model jShakura & Sunvaevll 19731: iNovikov & Thorn3 
Il973h which assumes that at a given radius, r, all dissipated en- 
ergy is released as radiation and that the emission at each ra- 
dius is locally given by a blackbody spectrum. The innermost 
stable circular orbit (ISCO), beyond which the particles in or- 
' bit become dynamically unstable and plunge into the black hole, 
. is taken to be a manifestation of the inner disk boundary. This 
' model, however, is only self-consistent in the limit at which 
the disk is geometrically razor thin. Consequently, its validity 
is restricted to luminosities below ~ 30% of the Eddington lu- 
minosity, Lndd = 1.26 X lO-'** (M/Mq) erg /s, when radiation 
pressure will cause the disk to be minimally inflated (see e.g. 
iMcClintock etani2006h . 

High luminosity, optically thick accretion disks are bet- 
ter described by " s lim" a ccretion disks models, introduced by 
lAbramowicz et al.] (Il988h and later elaborate d by several au- 
thor s, including the most re cent contribution of lS^dowskil (l2009t) 
and IS^dowski et al.l (l201 ih . who computed a network of fully 
relativistic slim-disk models (in Kerr geometry) that densely 
samples the relevant parameter space (accretion rate M, alpha 
viscosity a, Kerr spin parameter a,). Tables and routines to ex- 



tract information from the slim disk database are available on- 
linil 

LMC X-3 is a black hole (BH) binary system in the Large 
Magellanic Cloud (LMC) at a distance of 48.1 kpc (derived from 
Orosz et al.ii2009i) . It was discovered by the UHUR U satellite in 
1971 and reported to be a discrete X-ray source by iLeong et"all 
(Il97lh . The optica l counterpart of LMC X-3 was first identified 
as a faint OB star dWarren & Penfoldll 19751) and by subsequent 
spectroscopic studi es pinpointed as a B 3 V main-sequence star 
in a 1.7 day orbit dCowlev et al.lll983h. From the radial ve loc- 
ity shifts of H and He absorption lines ICowlev et alJ (Il983h de- 
rived a large mass function, /m - 2.3 + 0.3 Mq, which estab- 
lished LMC X-3 positively as a BH candidate with a BH mass 
7 Mq < Mbh ^ 14 Mq and the mass of the secondary star 
4 Mq < M2 < 8 Mq. They also pointed out the possibility that the 
donor star is not filling its Roche lobe. The source exhibits large 
intensity variations, mostly soft spectra and a low absorption 
column density along the line of sight. These properties make 
LMC X-3 an ideal laboratory for testing our understanding of 
accretion disk physics. Presently, the two binary parameters of 
interest to us are poorly known, namely the orbital inclination 
angle / and the mass M of the black hole (e.g., see Soria et alj 
^ 2001) . For the former, we adopt / = 66° (iKuip er e t al. \9m> 
and for the latter, the round value M - \QMq. Because the spin, 
fl,, derived via the continuum-fitting method depends strongly 
on these two uncertain parameters, the spin values we quote in 
this paper are highly uncertain and cannot be considered reliable 
estimates of the spin of LMC X-3 . 

In the continuum fitting technique, one determines the ra- 
dius of the inner edge of the accretion disc from the temperature 
maximum of the soft X-ray flux and assumes that this radius cor- 



http://users.camk.edu.pl/as/slimdisk 
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Fig. 1. Examples of local (measured in the disk co-rotating 
frame) flux profiles emitted from the surface of a slim accre- 
tion disk {solid lines) and Novikov-Thorne disk (dashed lines) 
models. Flux profiles are shown for three different mass accre- 
tion rates M = 0.03, 0.3, 0.6 M^m- Note the eff'ects of advection, 
clearly visible in the slim-disk curves, which start to be impor- 
tant at m > 0. 1 . 



responds to the innermost stable circular orbit of the black hole, 
from which then the black hole spin can be deduced (see th e 
pioneering works by IZhang et"al1ll997l; iGierlihski et al.ll2"001h . 
This method of measuring a, requires accurate estimates of the 
binary's black hole mass, distance and inclination. Given those, 
it depends only on the properties of the accretion disc model. 
The most recent studies of LMC X-3 by means of the contin- 
uum fitting method, based on Newtonian or relativistic thin disk 
models in the low luminosity limit, include measurements of the 
inner disk radius ( Steingr et al, 201( 3; Dunn et al. 201Q) and esti- 
mates of the spin ( Davis et al.l2006l) as wel l as an analysis of the 
intrinsic disk emission dKubota et al.ll2010l) . 

In this paper, we analyze hundreds of RXTE observations 
of LMC X-3 usi ng the newly developed spectral fitting rou- 
tine slimbb GBursa et al.ll201 ll) which is an improvement upon 
and successor to the commonly used models kerrbb/kerrbb2 
dLi et al.ll2005h to the extent that it allows to fit high luminos- 
ity data as it builds upon a more general, relativistic model of 
slim accretion disks and uses photon ray-tracing from the ac- 
tual location of the disk photosphere. We explore the properties 
of the binary system by analyzing fits for different viscosity pa- 
rameters. In Section |2] we introduce the main features of slim 
accretion disks and in Section|3]we describe the data set and the 
analysis. The results are discussed in Section|4] In Section|5]we 
present our conclusions. 
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Fig. 2. Calculated theoretical continuum spectra of both slim- 
disk {solid lines) and NT {dashed lines) models. The two panels 
show spectra calculated using either blackbody emission plus a 
constant hardening factor / = 1 .6 or using BHSPEC emission. 
Mass accretion rates for both models are slightly different from 
each other and they are chosen is such a way that both NT and 
SD spectra have the same integrated (observed) luminosity L - 
0.1, 0.3, 0.6LEdd- The observer is fixed at inclination 66°, the 
BH mass is 10 M© and the distance is 48.1 kpc. 



The location of the effective photosphere differs significantly 
from the equatorial plane with growing luminosity, although 
the relative vertical disk thickness is not large (H/r < 1), not 
even for the near-Eddington luminosities. Accordingly, our 
ray-tracing computations are taken from the actual disk p ho- 
tosphere to the observer at infinity (IS^dowski et al.ll2b09l) . 



2. Relativistic slim disks 

In the spectral fitting procedure we implement an accretion disk 
model based on the solutions of relativistic slim disk (SD) the- 
ory. The major improvement of this model over other com- 
monly used disk models is that it includes three effects that 
are dominant at high accr etion rates, and are not present in the 
iNovikov & Thornel (Il973h (NT) disk model: 

- Radial advection of heat, which plays a substantial role at 
higher luminosities, is present and significantly modifies the 
flux of radiation emitted at a given radius in the inner disk 
region. 

- The inner edge for the disk radiation deviates from ISCO at 
high luminosity and can be considerably closer to the black 
hole due to the advective transport of heat generated by vis- 
cous processes. 



In Figure [T] we present profiles of the emitted flux for a moder- 
ately rotating black hole (a, = 0.65) for three accretion rate^ 
m = M/MEdd = 0.03, 0.3 and 0.6, where Msdd = 2.23 x 
10'** M/Mq gs"' is the mass accretion rate corresponding to 
the Eddington luminosity of an accretion disk around a non- 
spinning black hole. It is clear that high mass accretion rates 
produce a significant rate of advection which results in a change 
of the emitted flux profile. In this regime, much of the heat gen- 
erated by viscosity is radiated from the portion of the disk which 
is closer to the black hole than in NT solutions. At the highest 
accretion rates, the radiative efficiency of the disk decreases, as 
some amount of thermal energy is never emitted because it is 
advected through the horizon. A detailed discuss ion of the prop- 
erties of the shm-disk models may be found in IS^dowski et al.l 
(1201 Ih . 
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2.1. Spectral model 

Based on the relativistic slim-disk model of IS^dowski et al.l 
( 201 ih. we have developed a new code, called slimbb 



Bursa et al.l 1201 lb. to be used with the X-ray spectral fitting 
tool XSPEC (lAmaudll996l) . Here we give only a brief summary 
of the model for completeness. For a chosen set of model pa- 
rameters (spin, mass accretion rate and c-viscosity), the slim 
disk solutions provide a number of physical quantities which 
are used by the slimbb code to construct integrated at-infinity 
spectra of the disk. These quantities are the radial profiles of 
the emitted flux from the photosphere, the surface density, the 
specific angular momentum, the fluid radial velocity and the lo- 
cation of the photosphere. Emission from the disk is ray-traced 
from the proper effective photosphere to an observer at infinity 
taking into a ccount all relativistic effects in the Kerr geometry 
(lBursall2006l) . and integrated over the disk surface to obtain the 
total observed continuum spectrum. While integrating the spec- 
tra, since the observed emission is affected by the opacity of the 
disk and can generally deviate from a pure blackbody, a color 
temperature correction (or hardening correction) factor needs to 
be applied to account for the non-blackbody contributions to the 
disk spectrum. This effect is is either incorporated using a disk 
atmosphere model (e.g., BHSPEC, iDavis & Hubenv 2006) or a 
single color correction factor is a ssumed for the effective tem- 
perature over the whole disk (e.g., iGierliriski & Donell2004[ and 
references therein). Our code offers both options, it can either 
integrate local blackbody spectra modified by a constant hard- 
ening factor (given as a parameter of the model) or it can di- 
rectly interpolate over l ocal disk annuli spectra constructed by 
iDavis & Hubenvl (l2006h . which already include all reprocessing 
effects on the emerging spectrum in the disk atmosphere, and 
integrate those to form the resulting spectrum without a need to 
specify an ad-hoc value of spectral hardening. Figure|2]compares 
"at-infinity" spectra of the slim disk and NT-disk models which 
are computed either using constant hardening or using BHSPEC. 
The two spectra of each NT/SD pair in the two panels are cho- 
sen such that they have the same luminosity and they all have the 
same fixed spin. In this comparison, clearly, the BHSPEC NT and 
SD spectra differ more from each other (especially at high en- 
ergies) than the same luminosity spectra with just a single value 
color temperature correction. 



3. Data reduction & analysis 

Our data have been obtained using the large-area Proportional 
Counter Array (PCA) aboard the Rossi X-ray Timing Explorer 
(RXTE). We use the same s ample of 712 PCU- 2 X-ray spec- 
tra as described in detail by ISteiner et al.l (1201 Ob . These spec- 
tra are mostly thermal and reduce in the course of the fitting 
procedure to 200-300 thermal spectra that obtain a sufficient 
goodness-of-fi{3(;t~/£/o/ < 2) and satisfy the following require- 
ments in the 2.55-25.0 keV energy band: (i) the photon index of 
the power-law component (see next paragraph) lies in a typical 
soft-state interval 2 < F < 3, (ii) the fraction of up-scattered 
photons is fuse < 10% and (iii) the ratio of disk flux to total 
flux is FdisklFtot ^ 60%. The analysis has been performed with 



the latest version of the X-ray spectral fitting software, XSPEC 
V. 12.6.1. 

Our full model combination in XSPEC reads TBabs * 
(simpl (g) slimbb): The photon absorption, TBabs, is deter- 
mined solely by the Galactic neutral hydrogen column density , 
«H, which we fix at the value 4 x 10^" cm"- (Page et al.,i2003l). 
The p ower-law component is generated by simpl (ISteiner et al.l 
I2009h . with a photon index, 1 < F < 10, and the fraction of up- 
scattered disk photons, fuse- The slim disk component, slimbb, 
is specified by ten parameters, two of which are free: we fit for 
the luminosity, 0.05 < lumin = Ldisk/i'Edd ^ 1-5, and the dimen- 
sionless BH spin, a. We freeze all binary system parameters (see 
Section[TJ: the black hole mass, M - 10 M©, the average distance 
to the LMC, D = 48.1 kpc, and the inclination, inc - 66°. We fix 
the viscosity parameter at either 0- = 0.1 or 0.01. At all times also 
fixed are the switch for limb darkening (on). If lag -I, the switch 
for vertical disk thicknes^ (on), vflag - 1, and since M, D and 
inc are known and fixed, the normalization is set to 1 . The switch 
for the spectral hardening factoid is set to /hard = -1, such that 
the code integrates local BHSPEC annuU spectra to produce the 
best total continuum spectrum. 



4. Results & discussion 

The standard thin disk model fails above roughly L a 0.3 Lndd 
due to the fact that a number of underpinning physical as- 
sumptions break down beyond this li mit. For this reason 
the spi n analysis with k e rrbb2, e.g . , in Shafee et al. ( 20061): 
McClint ock et alJ (l2006h : iGou et alJ (l2010h and ISteiner et al.1 
t2010i) has always been strictly limited to low luminosities and 
excludes thermal spectra with L > 0.3LEdd- The thin-disk model 
limitations become manifest when the mass accretion rate be- 
comes appreciably high, such that the disk is no longer able to 
radiate away all energy dissipated by viscous stresses. As a con- 
sequence the thin disk overheats and inflates and changes from 
being gas pressure dominated at low accretion rates to being ra- 
diation pressure dominated at higher mass accretion rates. The 
pressure perturbations set a viscous instability in motion and re- 
sult in the break down of the standard thin-disk model. An obser- 
vational consequence of the thin-disk model being used beyond 
its applicable limit (a; 30% Eddington; McClintock et al. 200^) 
could be that the spin appears to decrease with increasing lu- 
minosity. Naturally, the slim-disk model, which relaxes some of 
the most critical thin-disk assumptions, may provide a better de- 
scription of high-luminosity accretion disks and allow us to test 
the role of its free parameters. 

In the following, we compare the performance of our slim- 
disk model by fitting for the spin of the black hole in LMC X-3 . 
We compare the results with estimates obtained using the 
kerrbb2 model for two distinct values of viscosity. By mapping 
the respective spin-luminosity dependence we are able to assess 
(i) whether or not the slim-disk model performs better at high 
luminosities than the NT model and (ii) how viscosity affects 
the spin estimation. The influence of the a-viscosity parameter 
is then elaborated. 



^ Given the efficiency of conversion of rest mass to radiation en- 
ergy for spin a, = 0.65, these values approximately correspond to disk- 
integrated luminosities (oc J F(r)rdr) of Laisk= 0.1, 0.5 and 1.0 Z-Edd- 

^ The precise number of well-fitted thermal spectra found in the sam- 
ple is model dependent. 



The alternate setting, vflag = 0, starts the ray-tracing from the equa- 
torial plane. 

' Alternatively, /hard can be fixed at any value > 1 and the code inte- 
grates color-corrected local blackbody spectra. 
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Fig. 3. Comparison between the disk models slimbb (red diamonds) and kerrbb2 (orange triangles) for viscosity parameters 
a - 0.1 (left) and a = 0.01 (right). The dashed and dot-dashed lines represent linear fits to the data. In all cases we assume 
LMC X-3 contains a 10 M© black hole, is located at a distance of 48.1kpc and seen at an inclination of 66°. 



4.1. The spin of LMC X-3 

We apply the slim-disk model to our selected set of 
LMC X-3 spectra for two values of the alpha-viscosity parame- 
ter, a = 0.1 and 0.01. Then, we compare the results to analogous 
studies conducted with kerrbb2, the BHSPEC enhanced version 
of kerrbb, matching the slim-disk environment by switching off 
the returning radiation fla^ We present our results in Figure |3] 
in a form of comparative diagrams which map the spin against 
the luminosity. The error bars in all cases represent the 90% con- 
fidence interval. Error bars in luminosity which are smaller than 
the symbol width are suppressed. 

The left panel of Figure |3] shows that thin and slim disk re- 
sults with ff = 0.1 follow one another very closely and both 
models show an equally pronounced spin drop-off. Adopting the 
lower value, a - 0.01, in the right panel of Figure |3] delivers 
some improvement for slim disks. Although the situation is still 
not perfect, low alpha slim-disk models seem to be able to keep 
a constant spin over a larger span of luminosities, approximately 
up to 0.4LEdd- From this perspective, a lower value of viscosity, 
a<0.01, seems to be a more favorable choice. 

To compare the effect of different viscosity parameters and 
color correction treatment we then plot three different spin esti- 
mations based on slimbb. The bottom two sets in Figure|4]use 
the BHSPEC annuli spectra and show the effect of a, while the 
top set applies a constant hardening factor fixed at /hard = 1.6. 
This value is typical for the lower-luminosity RXTE spectra and 
is derived reversely from the best fit parameters of an arbitrary 
low-luminosity spectrum by freezing the best fit and thawing 
/hai-d instead. Note, that due to the fact that the spectral harden- 
ing is strongly degenerate with the BH spin one of them has to 
be fixed during the fitting procedure. 

Figure|4]shows that slimbb with a constant hardening factor 
produces the expected result: namely, that the same spin is ob- 



* Self-irradiation of the disk has in the case of a moderately rotating 
black hole only small impact on the spectrum, whi ch may be com pen- 
sated by adjusting the mass accretion rate in the fit (iLi et al.ll2005h . 
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Fig. 4. Three spin estimations with slimbb assuming a 10 M© 
black hole in LMC X-3 . Different alpha parameters and color 
correction treatments lead to a variety of trends with luminosity. 
The top two sets of fits have a low viscosity parameter, a = 0.01. 
The upper of the two (yellow circles) represents a fit using a 
constant hardening factor The middle pattern (orange triangles) 
shows a fit obtained using /hard calculated with BHSPEC. The bot- 
tom points (red diamonds) result from a high viscosity parame- 
ter, ff = 0.1, with BHSPEC hardening. 



tained over the full range of luminosity. We know, however, that 
the s pectral hardening factor generally increases with luminosity 
(e.g. lDavis & Hubenvll2006l) . Thus, holding /hard constant mim- 
ics the implicit luminosity dependence of another parameter 
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Fig. 5. Behavior of the o'-viscosity with luminosity for a fixed 
spin parameter value a, = 0.61. The shaded region represents a 
Gaussian fit to the 1-cr range of the viscosity parameter 

4.2. Viscosity and spectral liardening 

Viscosity aff'ects the spectra in such a way that for a given spin 
and luminosity a decreasing alpha value increases the local col- 
umn density, so that the flow is more optically thick and there 
is less effective up-scattering. Consequently, the emerging radi- 
ation has a softer spectrum. In other words, the lower the a, the 
lower the color correction of the spectrum. However, at a given 
viscosity, the hardening roughly scales as /hard L. Based on 
the fact that fits with a fixed hardening factor (although it does 
not represent a correct physical picture) can produce the desired 
result (i.e., that the resulting spin is independent of luminosity) 
we suggest that perhaps the alpha viscosity of the accretion flow 
changes with luminosity. In this sense, a slim accretion disk sce- 
nario with a variable viscosity parameter, which is anticorrelated 
with luminosity, could solve the spin drop issue. 

We can test this hypothesis with our slim-disk model by fix- 
ing all data at the mean spin value for the low-luminosity spec- 
tra a« 1/^0.2 LEdd - 0.61 and fitting the set of LMC X-3 spectra for 
0.001 < a < O.I. About 40% of the fits spread over the whole 
luminosity range pegged at the lower limit 0.001 and were ex- 
cluded. The remaining data points were fitted with a Gaussian to 
illustrate the preferred viscosity range. Figure |5]clearly indicates 
the trend towards ever lower viscosities as the luminosity in- 
creases. It also shows that in the low-luminosity limit, L < 0.25, 
quite in an agreement with theory, the viscosity parameter is de- 
generate and may assume any value. 

5. Conclusions 

5.1. kerrbb2 vs. slimbb 

In the previous Section we discussed the results of the contin- 
uum fitting method tested on the X-ray spectrum of the BH bi- 
nary system LMC X-3 . We used the disk models slimbb and 
kerrbb2 in combination with the power-law model simpl. For 
both disk models the qualitative behavior is similar Both mod- 
els are in excellent agreement in the low-luminosity limit. The 
slim disk solution, however, appears to be less strongly luminos- 
ity dependent than kerrbb2 (i.e., the drop in a, is less severe) - 
especially for low values of viscosity. Nonetheless, we conclude 
that applying the slim-disk model at a fixed value of a (which 
properly accounts for advection of heat at high luminosities) is 
insufficient to counter the high luminosity fall-off of the mea- 
sured spin. Some additional processes are likely to play a role. 
In the following subsections we discuss factors that may affect 
disk appearance and fitting results at high luminosities. 



5.2. How to produce softer spectra? 

Obtaining spin estimates consistent for all luminosities requires 
even softer theoretical spectra at high luminosity than obtained 
with the current slim-disk + BHSPEC model. Factors which 
so far may not have been taken into account "properly" in- 
clude: (i) global and luminosity dependent disk parameters, e.g., 
the viscosity parameter a, (ii) luminosity dependent disk atmo- 
sphere properties, i.e., the hardening factor /hard, (iii) wind mass 
loss l eading to the reduction of surface density (Takeuchi et al] 
I2OO9I) . 

Flux profiles of slim accretion disks for moderate and high 
accretion rates have more parameter dependencies than NT thin 
disks. While for the classical solutions of Novikov & Thome 
the disk emission is insensitive to the viscosity parameter a, 
this statement is no longer satisfied for slim disks at high ac- 
cretion rates. The form of the viscosity prescription, as well as 
assump tions about the disk ve rtical structure, influence the flux 
profile dS^dowski et al.ll20I ih . For the aP class of models the 
higher the value of a for a given mass accretion rate, the closer 
to the BH horizon the flux profile extends along the slim disk 
slope Teff oc r"'^^, reaching higher effective temperatures. The 
opposite limit, ff — > 0, defines the softest possible disk emission. 
Figure |4] shows that even with the lower value, a = 0.01, the re- 
quired softness at high luminosities cannot be obtained. Going to 
ever lower a values certainly helps, it is, however, not the proper 
way to tackle that problem. Other factors should also play a role. 

The BHSPEC routine accounts for photon up-scattering in the 
atmospheres based on radiative transfer calculations of plane- 
parallel disk annu li assuming a flat (oc p) dissipation profile 
dPavis et"ani2005h . Their results depend (through surface den- 
sity) on the assumed value of the a parameter - the higher the a, 
the harder is the emerging spectrum. Therefore, more consistent 
spin estimates may be obtained if one assumes that a changes 
with luminosity - if it dropped with increasing luminosity the 
outcome spectra would be softer 

Two other factors affecting Comptonization in disk atmo- 
spheres which are most profound at high luminosities have 
to be mentioned. Firstly, as recent radiative MHD simulations 
show, for high accretion rates significant wind outflow ap- 
pears to change the dynamical structure of the disk atmosphere 
which may affect the e ffectiveness of photon up-scattering 
(lOhsuga et alJ2005Ll2009h . However, such a mass loss also leads 
to reduction of the surface density, which is equivalent to in- 
creasing a, and thus may cause the opposite effect on the spec- 
tral hardness. Secondly, the higher the luminosity, the thicker the 
disk, and therefore, more radiation comes back as "returning ra- 
diation", which illuminates the upper disk atmosphere. This also 
may affect the Comptonized disk emission. 

5.3. The future of siim disks 

We saw in the previous sections that the new spectral model 
slimbb based on slim disks provides effective improvement 
and is a valid successor of the kerrbb/kerrbb2 models. 
Nonetheless, the general understanding of accretion disks is still 
far from being complete and the above mentioned issues are im- 
portant for all accretion disk models and need to be elaborated 
in future work. Considering slim disks in the high accretion rate 
regime is the first step towards a better understanding of accre- 
tion disks. 
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